Uncertainties in nuclear transition matrix elements for neutrinoless ββ\beta \beta decay II: the heavy Majorana neutrino mass mechanism

Employing four different parametrization of the pairing plus multipolar type of effective two-body interaction and three different parametrizations of Jastrow-type of short range correlations, the uncertainties in the nuclear transition matrix elements $M_{N}^{(0\nu)}$ due to the exchange of heavy Majorana neutrino for the $0^{+}\rightarrow 0^{+}$ transition of neutrinoless double beta decay of $^{94}$Zr, $^{96}$Zr, $^{98}$Mo, $^{100}$Mo, $^{104}$Ru, $^{110}$Pd, $^{128,130}$Te and $^{150}$Nd isotopes in the PHFB model are estimated to be around 25%. Excluding the nuclear transition matrix elements calculated with Miller-Spenser parametrization of Jastrow short range correlations, the uncertainties are found to be 10%-15% smaller